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ABSTRACT 

Massive quiescent galaxies at z > 1 have been found to have small physical sizes, hence to be superdense. 
Several mechanisms, including minor mergers, have been proposed for increasing galaxy sizes from high- to 
low-z. We search for superdense massive galaxies in the Wide-field Nearby Galaxy-cluster Survey (WINGS) 
of X-ray selected galaxy clusters at 0.04 <z< 0.07. We discover a significant population of superdense massive 
galaxies with masses and sizes comparable to those observed at high redshift. They approximately represent 
22% of all cluster galaxies more massive than 3xl0 10 M Q , are mostly SO galaxies, have a median effective 
radius (R e ) = 1.61 ±0.29kpc, a median Sersic index (n) = 3.0 ±0.6, and very old stellar populations with a 
median mass-weighted age of 12. 1 ± 1 .3 Gyr. We calculate a number density of 2.9xlO -2 Mpc~ 3 for superdense 
galaxies in local clusters, and a hard lower limit of 1.3xl0~ 5 Mpc~ 3 in the whole comoving volume between 
Z = 0.04 and z = 0.07. We find a relation between mass, effective radius and luminosity-weighted age in 
our cluster galaxies, which can mimic the claimed evolution of the radius with redshift, if not properly taken 
into account. We compare our data with spectroscopic high-z surveys and find that - when stellar masses 
are considered - there is consistency with the local WINGS galaxy sizes out to z ~ 2, while a discrepancy 
of a factor of 3 exists with the only spectroscopic z > 2 study. In contrast, there is strong evidence for a large 
evolution in radius for the most massive galaxies with M* > 4x10' 'M compared to similarly massive galaxies 
in WINGS, i.e. the BCGs. 

Subject headings: galaxies: clusters: general — galaxies: evolution — galaxies: structure — galaxies: funda- 
mental parameters 
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1. INTRODUCTION 

In the last years high-z studies have uncovered a 
considerable number of massive galaxi es with relative! 
small effective rad i i (see , among others. iDaddi et al.1 1 20051 
Truiillo et all [20061 [2007t iToft et all 12007 1 : IZirm et all 1200' 
Buitra goet all 120081: ICimatti et all l2008t Ivan Dokkum et alj 
20081: ISaracco et all l2009t Ivan der Wei " et all 120081) . or. in 
other words, superdense galaxies 1 (hereafter, SDGs). Al- 
though the datasets and methodologies are quite different, 
they all agree on the fact that a population of massive and 
compact passive galaxies at z > 1 does exist, with sizes a fac- 
tor of at least 3 less that their low-z counterparts of the same 
mass. 

Other studies, using different samples extracted from 
SDSS, have found a complete absence of such galaxies 
with old stellar pop ulation ages (see, e.g.. IShen et al.l [2O03h 
Truiillo et al. 2009). This implies the necessity of an evolu- 
tion in radius with redshift, and it is often considered a proof 
that these high-z SDGs have undergone significant (minor 
and/or major) me rging events along their histories. 

More recently, iTavlor et alj (120091) identify 63 M* > 5 x 
10 10 M Q red sequence z ~ 0.1 SDG local candidates which 
are smaller than the median mass-size relation by a factor 
of at least 2. These local galaxies have sizes which are 
compatible with many of the high-z ones, but they seem 

As far as physical densities are concerned, the galaxies under investi- 
gation here are thou ght not to be extreme (see, e.g.,|Bezanson et al] |2009t 
Hopkins et al. 20091) 



2) than the most distant 
m galaxies recently discov- 



anyway larger (by a factor of ' 
z > 2 and massive M* > lO j'M 

ered (see, e.g. iBuitrago et alj|2008t Ivan Dokkum et all 120081: 
iDamianov et al.ll2009l) . 

Different scenarios have been proposed to explain the com- 
pactness of these galaxies and their subsequent evolution. 
One of the simplest ideas is to assume that high-z masses 
and/or R e measurements, are incorrect. It is true that sys- 
temati c effects can easily poll u te these measurement s (see, 
e.g., Ivan Dokkum et al] 12008: Be zanson et al.l 12009). such 
as low signal-to-noise ratios, limitations in resolution, uncer- 
tai nty on the IMF, models an d SED fitting, etc. A recent study 
of Ivan Dokkum et"ai1 d2009l) . analyzing a 29 hours exposure 
spectra of one of their 9 high-z SDGs, has found a very high 
velocity dispersion (^500km/s), which is consistent with its 
compact nature and stellar mass from SED fitting. Obviously, 
this is only one case, and it must be confirmed, but it gives an 
indication that measurements errors may not be the explana- 
tion for the existence of these high-z mas sive compact objects. 

On the other hand, Ma ncini et al.l d2009l) . over the 2 degrees 
COSMOS field, select 12 quiescent massive galaxies at z ~ 
1.5, and found sizes mostly compatible with the local mass- 
radius relation. Based on mocked images of high-z galaxies 
initially laying on the local mass-radius relation, they claim 
that size measures performed on low S/N images are likely to 
give systematically lower R e . 

By studying a sample of Brightest Cluster Galaxies 
(BCGs), based on considerations on the evolution of the mass- 
function with redshift and the merging rates in numerical sim- 
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ulations. iBernardil (12009) claims that the main mechanism to 
let galaxies increase their radius without gaining too much in 
mass, is through minor mergers. This is in general interpreted 
as an independent indication that very old galaxies have in- 
creased their size during their evolution by means of minor 
merger. However, these findings are not conclusive, as it is 
not clear whether mergers alone can efficiently puff up galax- 
ies by the required amount, as many parameters are involved 
in this mechanism (see, among othersjK hochfar & Silkll2006t 
Ivan der Wei et al.ll2009t Houng et al.ll2009l) . Furthermore, this 
mechanism could be efficient only for BCGs, becoming less 
and less relevant for other kind of galaxies. 

Another viable explanation for puffing up galaxies might 
be connected with the quasar phase that the se galaxies hav e 
likely undergone from z ~ 2.5 to z ~ 1 (Fan et al.l [2008b . 
which caused a dramatic mass loss with consequent expan- 
sion of the galaxy. However, such expansion would take place 
at most in a few dynamical times (~ 8xl0 8 Gyr), causing only 
a few systems to be caught as quiescent and still compact (see 
iMancini et all2009l) . 

Whatever the evolution mechanism, local clusters could be 
an ideal place for SDGs, as they probably reside in very dense 
environments at high redshift too; this is supported by the 
strong clusterin g (R a 8 — lOMpc ) of quiescent 2 compact 
high-z galaxies (Ci matti et alj|2008l and references therein). 
At least a fraction of these objects may have survived till 
recent cosmic epochs; some models predict that 10% of 
galaxies have had no significant transformations since z ~ 2 
(Hopk ins et al.ll2009bl and references therein). Hence, it is 
plausible that a certain number of very old compact galaxies 
are found in local galaxy clusters. I n this contest the W lde- 
field Nearby Galaxy-clusters Survey (Fasano et al. 2006|) is a 
suitable survey where to search for such objects. 

The layout of this paper is as follows. In §2 we describe the 
data set we used to search for SDGs. In this section particular 
care is given to the homogenisation of our data with litera- 
ture data, and some caveats are discussed. In §3 we present 
our sample of compact galaxies. In §4 we discuss the com- 
parison with high-z data and the selection effects that occur 
when considering old stellar populations at high-z, and give 
number densities and frequencies. In §5 we discuss the lo- 
cal SDSS mass-radius relation used by high-z studies as a lo- 
cal reference and some reasons why recent works may have 
missed local counterparts to massive high-z compact galax- 
ies. In §6 we describe in more detail all photometric, spec- 
troscopic and intrinsic properties of our compact sample, and 
finally we draw our conclusions in §7. 

Throughout this paper we will use the cosmology (Hq, Q, m , 
A ) = (70,0.3,0.7). 

2. THE DATA SET 

The galaxies examined in this paper are part of the Wlde - 
field Nearby Galaxy-clusters Survey (Fasano e t al.l 120061) . 
WINGS 3 is a multiwavelength survey especially designed to 
provide the first robust characterization of the photometric 
and spectroscopic properties of galaxies in nearby clusters, 
and to determine the variations of these properties as a func- 
tion of galaxy mass and environment. 

Clusters were selected in the X-ray from the ROSAT 

2 High-z galaxies are considered quiescent when their luminosity- 
weighted age is > 1.5Gyr; we call quiescent WINGS cluster members with 
luminosity-weighted ages > lOGyr, i.e. quiescent atj ~ 1.5 

3 Please refer to WINGS Website for updated details on the survey and its 
products, http : / / web . oapd . inaf . it/ wings 



Brightest Cluster Sample and its extension (Ebeling et al. 
1998, 2000) and the X-ray Brightest Abell-type Cluster sam- 
ple (Ebeling et al. 1996). WINGS clusters cover a wide 
range of velocity dispersion cr c ius> typically between 500 and 
1100 kms~', and X-ray luminosity Lx, typically 0.2 — 5 x 
10 44 erg/s. 

The survey core, based on optical B, V imaging of 78 n earby 
(0.04 < z < 0.07) galaxy-clusters dVarela et alj 120091) . has 
been complemented by several ancillary projects: (i) a spec- 
troscopic follow up of a subsample of 48 clusters, obtained 
with the spectrog raphs WYFFOS@WHT and 2dF@AAT 
(ICava et al.l 120091) : (ii) near-infrared (J, K) imaging of a 
subsample of 28 cluster s obtained with WFC AM @ UKIRT 
dValentinuzzi et al]20 09) ; (iii) U broad- and H Q narrow-band 
imaging of subsamples of WINGS clusters, obtained with 
wi de-field cameras at di fferent telescopes (INT, LBT, Bok, 
see[Omizzolo et al. 2009, in preparation). 

In the following, we will use only spectroscopically con- 
firmed members of the subset of WINGS clusters that have an 
average spectroscopic completeness larger than 50% (21 out 
of 78 clusters). Our completeness is essentially independent 
of distance to the cluster center for most clusters, and is com- 
pletely independent of galaxy radius. The only cri terion used 
for sp ectroscopic selection was galaxy magnitude dCava et ail 
2009), but given that separate configurations were used to take 
spectra of bright and faint galaxies and also due to fiber colli- 
sion effects, completeness turns out to be rather flat even with 
magnitude for most clusters. 

In this paper WINGS results are compared with literature 
data at 0.9 < z < 2.5. Several studies have investigated the 
sizes of distant quiescent galaxies, but we only consider here 
high-z datasets based on spectroscopic redshifts that give 
high quality masses and s izes, while other works that used 
photometric redshift s (i.e.. iToft et aTll2007t IZirm et alJl2007t 
iBuitrago et al.l 120081) are not included i n the present study . 
We use th e following dataset s: HUDF dDaddi et al.l l2005b 
MUNI CS dTruiillo et al] 120 06). MUSY C (Ivan Dokkum et all 
1 20081). ISaracco et al] d20 09t). GMASS (ICimatti et al.l 12008). 
Ivan der Wei et al] (120081) and iDamianov et al] d2009l) . The 
data, methods of analysis and, most importantly, selection cri- 
teria for these samples clearly differ from one study to an- 
other. In the comparison amongst different samples it is of 
paramount importance to account for differences in models 
and IMF adopted (see following se ctions). We stress t hat all 
of these works, with the exception of Ivan der Wei et aT] (120081) 
that have used a visual early-type morphological classifica- 
tion, have selected their galaxies to have already old (typi- 
cally 1.5-2 Gyr) stellar populations at that redshift based on 
SED spectral fitting, line index age dating, absence of sig- 
nificant emission lines, or other spectro-photometric analysis 
methods. 

2.1. Surface Photometry and Morphology in WINGS 

WINGS effective-radii, axial ratios and Sersic indexes 
are measured on the V-band images wi th GASPHOT 
dPignatelli et al] 120061; iD'Onofrio et al]|2009l) . an automated 
tool which performs a simultaneous fit of the major and mi- 
nor axis light growth curves with a 2D flattened Sersic-law, 
convolved by the appropriate, space-varying PSF In this way 
GASPHOT exploits the robustness of the ID fitting technique, 
keeping at the same time the capability (typical of the 2D ap- 
proach) of dealing with PSF convolution in the innermost re- 
gions. 

GASPHOT has proved to be very robust in recovering the 
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best fitting parameters, and to give the appropriate weight to 
the external parts of the galaxies, where PSF effects are neg- 
ligible. Indeed we tested GASPHOT on more than 15,000 
simulated and real galaxies, obtaining robust upper limits for 
the errors of the global parameters of galaxies, even for n on- 
Sersic profiles and blended objects (Pig natelli et al]l2006l) . 

GASP HOT was also tes ted against the widely used tools 
GALF IT dPeng et al.ll2002h an d GIM2D (Marlea u^Smwdl 
1998): it has been found (see. iPignatelli et al.l 120061 section 
6) that the performances of these tools are quite similar for 
large and regular simulated galaxies, while GASPHOT has 
proved to be more robust for real galaxies with some kind of 
irregularity or blending, which is a crucial feature when deal- 
ing with blind surface photometry of huge galaxy samples. In 
section 13.11 we show a comparison of GASPHOT estimates 
with literature data. 

The GASPHOT output effective radius R e value is calcu- 
lated along the major-axis, and for the purposes of this paper 
is circularized with the usual formula: 

R( drC ) = R^j™) ■ yjtyl (1) 

where a and b are the major- and minor-axis of the best-fit 
model, respectively. 

WINGS morphologies are derived fro m V images using 
the purposely devised tool MORPHOT dFasano et al.ll2009l 
in preparation). Our approach is a generalization of the non- 
para metric method p roposed by Conse lice et al.l (120001) (see 
also, IConse lice 2003). In particular, we have extended the 
classical CAS (Concentration/ Asymmetry/clumpinesS) pa- 
rameter set by introducing a number of additional, suitably 
devised morphological indicators, using a final set of 10 pa- 
rameters. A control sample of 1,000 visually classified galax- 
ies has been used to calibrate the whole set of morphological 
indicators, with the aim of identifying the best sub-set among 
them, as well as of analyzing how they depend on galaxy 
size, flattening and S/N ratio. The morphological indicators 
have been combined with two independent methods, a Maxi- 
mum Likelihood analysis and a Neural Network trained on the 
control sample of visually classified galaxies. The final, au- 
tomatic morphological classification combines the results of 
both methods. We have verified that our automatic morpho- 
logical classification reproduces quite well the visual classi- 
fication by two of us (AD and GF). In particular, the robust- 
ness and reliability of the MORPHOT results turn out to be 
comparable with the typical values obtained comparing each 
other the visual classifica tions obtained by d ifferent (experi- 
enced) human classifiers (Fas ano et al.l l2009, in preparation). 
Although MORPHOT provi des a fine classificatio n following 
the "Revised Hubble Type" de Vaucoul eursl (1 19741) . we will 
use in the following just three broad morphological classes, 
ellipticals, SOs (together early-type) and late-type, where the 
late-type class includes any galaxy later than an SO. 

2.2. Stellar masses, ages and metallicity 

Stellar masses of WINGS galaxies have been determined 
by fitting the optical spectrum (in the range ~ 3600-^ ~ 7000 
A), with the sp ectro-photometric model fully described in 
iFritz et all (12007b . All the main spectro-photometric features 
(such as the continuum flux and shape, and the equivalent 
widths of emission and absorption lines) are reproduced by 
summing the theoretical spectra of Simple Stellar Population 
(SSP) of 13 different ages (from 3xl0 6 to - 14xl0 9 years). 

Dust extinction is allowed to vary as a function of SSP age, 
in a screen uniformly distributed in front of the stars. The 



Galactic extinction law follows ICardelli et al.l (1 19891) scheme , 
with R v = 3.1. As explained in detail in IFritz et all d2007l) . 
for the fit we use a fixed metallicity, exploring three values: 
Z=0.004, Z=0.02 and Z=0.05. The adopted star formation 
histories and stellar masses refer to the model with the metal- 
licity value that provides the lowest x 2 - The lowest x 2 f° r the 
great majority of our super-dense galaxies yields either solar 
or supersolar metallicities. 

SSP spectra are built using Padova evolutionary tracks 
and the observed MILES spect r al libr ary (ISanchez-BlazquezI 
2004: ISanchez-Blazquez et alj 120061) for ages older than 
10 9 vr. complemented by the lJacobveT al. (1984) library for 
young SSPs, and in the UV and infrared by means of the Ku- 
rucz theoretical library. Nebular emission is also included, 
modeled with values that are typical of Hll regions: this sig- 
nificantly affects spectra of SSPs younger than ~ 2xl0 7 years. 

2.2.1. Stellar masses 

There are three main definitions of galaxy stellar 
mass derived by means of spectral synthesis (see, e.g., 
lLonghetti & Saraccoll2009l for details): 

1. the initial mass of the SSP, i.e. the mass of all the SSP 
stars at the moment of their formation. This mass does 
not depend on the SSP's age, being fixed once and for 
all; 

2. the mass locked into stars, including stellar remnants, 
at any time; 

3. the mass of stars that are still in the nuclear burning 
phase {i.e. no remnants included), at any time. 

The difference between these definitions (up to a factor of 2 
from definition 1 to definition 3, depending on several model 
parameters, such as e.g. the IMF) is a function of the SSP 
age, as the fraction of gas which is returned to the interstellar 
medium and the fraction of stars that evolve into remnants 
increase with time. In this paper we only use masses derived 
from definition 2. 

Our spectra are taken within a 2 arcseconds aperture fiber. 
For the purpose of computing total stellar masses and star for- 
mation rates, model spectra are rescaled to match the observed 
total V mag nitude, i.e. the SExtractor MAG_AUTO from 
IVarelaetalJ (l2009). This assumes that color gradients within 
galaxies are negligible. In order to take into account color 
gradient effects, we apply a correction of A(B — V) dex 4 to 
o ur masses. This correct ion is based on the prescription given 
in lBell& de Jon3d200l . 

2.2.2. Mass- and luminosity- weighted ages 

From our spectral analysis, it is possible to derive an es- 
timate of the avera ge age of the stars in a gal axy. Follow- 
ing the definition of lCid Fernandes et alj (120051) . we compute 
the luminosity-weighted age by weighting the age of each 
SSP composing the integrated spectrum with its bolometric 
flux. This provides an estimate of the average age of the 
stars weighted by the light we actually observe. A mass- 
weighted age is computed in a similar way: each SSP age 
is weighted with its mass value. The mass-weighted age is 
the "true" average age of the galaxy's stars. For our sample of 

4 A(B - V) = (B - V) flb er — (B- V)iokpc, where lOkpc is the physical 
aperture diameter. The median correction in mass for galaxies with M* > 
3xl0 10 M© is -0.05dex 
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cluster galaxies mass-weighted ages are systematically larger 
(~ 2Gyr) than luminosity-weighted ages. 

2.2.3. IMF and model differences 

The spectro-photometric analysis performed on the 
WINGS spectra was done assuming a Salpeter (1955) IMF 
with masses in the range 0.15 -f- 120 M Q . We then rescale 
both o ur values and all those from the literature to a iKroupal 
(2001) IMF, with masses in the range 0.01 50 M . It is also 
extremely important to properly match models that use dif- 
ferent treatments of the thermally pulsating asymptotic giant 
branch phase (TP-AGB). When needed, i.e. for high-z liter- 
ature galaxies, all the mass va lues were r escale d in order to 
match those obtained with the Ma raston] (12005b models, ap- 
plying a correction of 0.15 dex to the masses derived from 
Bruzu al and Chariot models, as prescribed by Cima tti et alj 
(2008). This difference in mass is strongly depending on 
the stellar population age, becoming practically negligible for 
ages older than 2/3Gyrs, i.e. it is ininfluent for local cluster 
galaxies, which are extremely old. 

3. RESULTS: WINGS SUPERDENSE GALAXIES 

In the bottom panel of FigQ] we plot the circularized effec- 
tive radius R e as a function of stellar mass for spectroscopi- 
cally confirmed WINGS cluster members with stellar masses 
M» > 10 10 Mq. In the upper-panel, we plot the mean mass 
surface density inside R e : 

as usually defined by other authors (see , among others, 
ICimatti et al.ll2008tlvan der Wei et al.ll2008b . 

After rescaling all the masses to the same IKroupal d200Th 
IMF and models (see previous section), we overplot high- 
z data from the literature as large open symbols. It can be 
seen that literature data cover a large range of masses and 
radii/densities, and a considerable fraction doesn't even reside 
in the highest mass-density locus. 

We choose the SDGs WINGS subsample (larger blue stars 
and red dots) in order to match as much as possible the po- 
sition of high-z SDG data in FigQ] (region inside the dashed 
lines), applying the following criteria: 

3xl0 10 M Q <M* <4xlO n M (3) 
E 5 o>3xl0 9 M Q kpc- 2 (4) 

We exclude from the SDG sample the Brightest-Cluster- 
Galaxies and other galaxies more massive than 4x10' 'M Q , 
as they may have a more complex formation history (dry /wet 
minor/major merger) which could, in principle, pollute our 
analysis (see §5 and Fig|5j. 

The resulting sample consists of 134 galaxies. Analyzing 
individually both images and spectra of this sample, we de- 
cided to exclude 12 of them because of close companions, 
bad chip regions, or low S/N spectra. From here on we re- 
fer to the remaining 122 objects as the WINGS SDGs sam- 
ple, that includes nearly 22% of all cluster members with 
M* > 3xl0 10 M Q ; we recall here that we are using only the 
spectroscopic confirmed members of a subset of 21 (all from 
the southern emisphere) out of 78 WINGS clusters. Regard- 
ing their morphologies, 31 of them are ellipticals, 78 SOs and 
13 late-type galaxies. 

In FiglJwe also draw the lShen et all d2003l) SDSS low-z re- 
lation for early-type galaxies, selected to have a Sersic index 



n>2.5 (full black line, with lcr and 2a as dotted lines), com- 
monly used by high-z studies as a reference point for the local 
mass-radius relation. We note at this point that the WINGS 
SDGs sample is found at more than 2a from the mean SDSS 
sample. In §5 we will discuss the local SDSS relation in more 
detail. 

It is clear that in the WINGS dataset we do find a consid- 
erable number of galaxies with masses, radii and mass den- 
sities typical of high-z SDGs. The only high-z samples that 
stands out for their extreme densities and lo w radii are 6 of 
the 9 galaxies from Ivan Dokkum et al. (2008) and 2 of the 10 
galaxies from iDamjanov et al.1 (120091) that do not have local 
WINGS counterparts. We will discuss further these cases in 
the following section. 

We note that for masses M* < 3xl0 I0 M Q there is a sig- 
nificant decrease in the frequency of SDGs (FigsQ] and 0: 
we checked whether this is due to completeness effects 
that could result in systematically missing effective radius 
and/or morphology measurements of small objects for low 
masses/luminosities, but this is not the case. We speculate 
that this rapid decrease in number is an indication that a mini- 
mum threshold in mass is required to form compact galaxies. 

3.1. Comparison of sizes and masses with literature data 

It is widely known that stellar masses based on spectro- 
photometric models have typical errors of ~ 0.2 dex. A 
crucial issue is to ensure that low-z masses are comparable 
to high-z ones. Dynamical masses from integral field spec- 
troscopy or virial masses from central velocity dispersions 
can, in principle, be the so lution for this kind of studies (see, 
e.g. Ivan der Wei et al.l 120081) . but it is very time consuming 
and very difficult to apply at z > 2. 

We calculated the virial masses of our WINGS galaxies 
from the velocity dispersions we found from literature: they 
are ^0.14 dex heavier than our Kroupa '01 IMF stellar 
masse s (the same offset was recovered by Carjrjellari et al. 
2006). Some 20% of our galaxies (either compact or nor- 
mal ones) present an excess of stellar mass. This can be ex- 
plained by the large uncertainties involved and by the need 
of accurate dyn a mical models, as thoroughly explained in 
ICappellari et al.l (120061) . Integral field spettroscopy of our 
SDGs will further clarify this issue. 

The SDGs central velocity dispersions a„ range from 100 
up to SOOkms- 1 ((<r ) = 180 ± SOkms" 1 ). These values 
are signific antly smaller than t he val ue of ~ 500kms _1 pre- 
sented by van Dokk um et alj d2009l) for a compact galaxy 
at z = 2.2 drawn from their MUSYC sample. Instead, 
they are s urprinsingly in agreement with measurements at 
z> 1-5 bylC enarro & Trujillo (2009), recently confirmed by 
ICappellari et al.1 d2009l) on GMASS galaxies in the redshift 
range 1.4 <z< 2.0. 

In the the top panel of Fig|2]we show the comparison be- 
tween WINGS early-type effective radii estimates and both 
NYU-VAGC 5 (Sersic fit, left panel, based on the spectro- 
scopic SDSS-DR7 catalogs) and SDSS 6 (de Vaucoleurs fit, 
right panel, based on the photometric SDSS-DR7 catalogs) 
values. We plot only early-type WINGS spectroscopic mem- 
bers with M* > 3xl0 I0 M Q , and for this comparison we use all 
WINGS survey data (not only the 21 clusters of this paper). 
Big coloured dots are the WINGS SDGs in common, red for 
ellipticals and green for SO's. It is apparent that most of the 

5 http://sdss.physics.nyu.edu/vagc/ 

6 http://cas.sdss.org/dr7/en/tools/crossid/crossid.asp 
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FIG. 1 . — The circularized effective radius R e and the mass-density inside R e as a function of stellar mass for all WINGS galaxies with M* > 1O 1O M , for 
the subsample of 21 clusters considered (see text). Blue and red tiny dots are late- (later than SO) and early-type (ellipticals and SOs) WINGS cluster galaxies, 
respectively. The region corresponding to our SDGs definition is delimited by the dashed lines in the top panel. The corresponding larger blue stars and red dots 
mark the WINGS SDGs. The black solid line is the SDSS-DR4 Shen et al. (2003) relation with dotted icr and 2cr lines. Open symbols are SDGs from high-z 
studies, see text for references. 



points lie in the region within 0.1 dex difference, and the me- 
dians do not show important offsets. There are some (~ 5%) 
that seem to show an underestimation of GASPHOT radii 
with respect to NYU-VAGC and SDSS. We checked "by eye" 
all these cases and found that many of these present a disturb- 
ing star in proximity of the galaxy, have a close galaxy pairs 
or groups, or a galaxy size below (at the limit of) the SDSS 
resolution. Due to the robustness of GASPHOT in these pe- 
culiar conditions, and to the deep er and higher- resolution im- 
ages of WINGS (for details, see iFasano et al.l l2006). we are 
tempted to think that our estimates are more reliable even in 
those cases. More interestingly, SDGs galaxies have com- 
patible R e measurements, as only one case shows a WINGS 
radius much lower than the literature one. Again, this is an SO 
strongly contaminated by a secondary object. 

In the bottom panels of Figj2] we show that our masses 
are in good a g reeme nt with the SDSS-DR4 estimates from 
iGallazzi et al.l ( 12005b . even though the scatter is high (~ 
0.15 dex), while we find an offset of - 0.09 dex with SDSS- 



DR7 masses 7 . Here we consider only the 21 clusters used in 
this paper (all from the southern emisphere), that have both 
high quality photometry and reliable masses based on high 
quality spectroscopic data. Only 3 of the 21 clusters are in 
common with SDSS (namely, A2399, and, partly, A119 and 
A957x which are covered at the 40% level), therefore the 
number of SDG galaxies that can be used for this comparison 
is small. Our mass estimates agree with masses calculated 
with (B-V) col o r and total V band fluxes with the recipe of 
iBell & de Jongl (|2001|). We note h e re tha t IShen et al l (120031) 
uses SDSS-DR4 iKauffmann et all (2003) masses which are 
lo wer by a factor of 0.0 7dex and of 0.13dex when compared 
to IGallazzi et al.l d2005l) and our masses, respectively. 

As a furthe r check we have run t he popular Hyperzmass 
software (see iBolzonella et"afl 120001) on B,V filters, finding 
no appreciable offset with WINGS masses. 

More interestingly, these comparisons and tests show that 
the SDGs are not systematically extremes, giving more reli- 
ability to our conclusions on compactness of massive local 



7 http://www.mpa-garching.mpg.de/SDSS/DR7/Data/stellarmass.html 
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FIG. 2. — Consistency check for WINGS sizes and masses. Only early- 
type WINGS spectroscopically confirmed member galaxies are shown. Big 
coloured dots are SDGs, red for ellipticals and green for SO's. Top panels: 
difference in R c estimates (dex) between GASPHOT, NYU-VAGC (Sersic 
fit, values from spectroscopic DR7 catalogs) and SDSS-DR7 (De Vaucoleurs 
fit, from photometric DR7 catalogs) values in arcseconds (see text for de- 
tail s). .Bottom panels: d ifference in mass (dex) between WIN GS estimates 
and lGallazzi etall f2005), SDSS-DR7 and Bell & de Jond<200i1) . 

WINGS cluster galaxies. 

4. COMPARISON WITH HIGH-Z: THE IMPORTANCE OF STELLAR 
AGE SELECTION EFFECTS 

Though, as shown in the previous section, there are low-z 
superdense counterparts to the high-z SDGs, and they repre- 
sent a significant fraction of the local cluster massive galaxy 
population, it remains to be addressed whether the prevalence 
of SDGs among the observed high-z galaxies requires a size 
evolution in a significant fraction of all massive galaxies. 

In Figf3]we show the combined effect of stellar mass and 
effective radius in determining the stellar age of a galaxy. 
The grey scale corresponds to luminosity-weighted ages (see 
legend for details): at fixed mass, smaller galaxies are older, 
while for a fixed radius, more massive galaxies are older. The 
same general trend is preserved if mass-weighted ages are 
used, so that even the formation epoch of the bulk of the stel- 
lar mass of these objects depends simultaneously on stellar 
mass and radius. 

Compari ng the sizes of massive high-z galaxies with 
the SDSS She n et all (120031) relation, several authors have 
claimed the necessity of an evolution of the size of such galax- 
ies with redshift, at least of a factor of 3 (0.5dex). While 
we will address the necessity to "properly" calibrate the lo- 
cal mass-radius relation in the next section, we want now to 
focus on the effect of a luminosity-weighted age selection. 
All high-z studies shown in Fig.[T]have selected their galaxies 
to be "old" on the basis of their stellar population properties 
(either SED fitting, lack of significant emission lines, spectral 
features etc), which translates into selecting galaxies with a 
luminosity-weighted age at least 1.5-2 Gyr old at the redshift 
they are observed. 

In Fig|4]we show the median R e of WINGS galaxies (filled 



dots) with luminosity- weighted ages older (by > 1 .5Gyr) than 
the age corresponding to the plotted redshift. The three panels 
refer to three stellar mass intervals, chosen to match the vari- 
ous high-z samples and to have a sufficient number of galaxies 
for statistics. The mean sizes of high-z literature data are plot- 
ted in color symbols (see legend). The dotted points have less 
than 3 galaxies, while all other have at least 3 galaxies in each 
interva l. The magenta pentagons are the Ivan der Wei et al.l 
(2008) data, who use virial masses on the basis of central 
velocity dispersion measurements; we correct these masses 
with a mean 0. 15dex contribution of dark matter in order to be 
globally compa tible with our ste l lar ma sses . Even with this 
correction, the Ivan der Wei etail (12008) "stellar-like" masses 
we derive may not be completely consistent with the stellar 
mass estimates of all other samples used i n this paper, so the 
comparison with Ivan der Wei et al.1 (120081) should be treated 
with caution. 

On top of the main well-known correlation of radius with 
mass (more massive galaxies are on average bigger in size), 
we find a noticeable decrease of the median radius with in- 
creasing redshift when galaxies are selected to be old at that 
redshift. The older the stellar population is selected, the 
smaller the median effective radius. 

Stellar age selection effects are therefore important: high-z 
studies find preferably compact galaxies because they select 
them to be old (i.e., to have only old stars). Assuming the 
sizes of today cluster massive galaxies to be representative of 
the sizes of all massive galaxies regardless of environment, 
we speculate that if high-z studies would include galaxies 
of all luminosity-weighted ages (young and old), they would 
find median effective radius values compatible with the global 
WINGS mass-radius relation, at least as far as z ~ 2 (i.e. our 
limit of resolution in age). 

In other words, a galaxy which is seen to be sta r-forming 
(and with larger sizes, see, e.g., iKriek et al.ll2009l) at, say, 
Z ~ 2 will obviously not be included in the high-z passive 
galaxy samples, which will have a smaller size. As we reach 
z ~ 0, those larger star-forming high-z galaxies that have be- 
come passive in the meantime, cause the local median size 
of passive galaxies to be apparently larger than that at high-z. 
The correct thing to do is to compare those high-z sizes with 
the sizes of local galaxies which were quiescent at those high 
redshifts. 

Fig|4] shows that, in general, no evolution in radius is re- 
quired for most of the high-z samples we consider in the 
present study. The majority of high-z datapoints are consistent 
within la with the WINGS estimates, even when the smallest 
errorb ars (errors on the means) are considered. 

The Ivan der Wei et al.l (|2008) datapoints, which we recall 
are obtained from dynamical masses, and therefore on stel- 
lar mass estimates that are not homogeneous with all other 
samples, show at most a factor 1.2 to 1.5 of evolution in size, 
much lower than the factor 3 claimed in th e literature. 

We note that when lSaracco et all d2009) divide their high-z 
sample in two classes of galaxies characterized by old (Lw- 
age ~ 3.5Gyr) and young (Lw-age ~ 1.5Gyr) stellar popula- 
tions, they find younger galaxies to have sizes compatible with 
the local mass-radius relation, in contrast with older galaxies 
that have smaller radii. The dependence of galaxy stellar age 
on galaxy size, at a fixed mass, must clearly be already estab- 
lished at z ~ 1.5. Our results confir m recent findings of other 
authors: Shank ar" fe BernardH d2009) have recently pointed out 
that, at fixed mass, smaller galaxies have older s tellar popu- 
lations, and the same conclusion is implicit in iGraves et al.l 
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FIG. 3. — Same as FigfT]but only for WINGS galaxies and with grey-scale filled dots showing luminosity-weighted ages. The black color is assigned to 
luminosity-weighted age > 9Gyr (equivalent to beeing quiescent at z ~ 1.3). More massive galaxies tend to have older ages and, at a given mass, galaxies with 
smaller radii are older (see §4). 

We tested the robustness of our conclusions by using 
masses with BC03 models instead of MA05, and luminosity- 
weighted ages calculated for a fixed metallicity (either solar or 
supersolar) for all galaxies. In all these cases the conclusions 
drawn from Fig|4]remain the same. 

Our conclusion is chall enged by the lowest median sizes of 
Ivan Dokkum et al] ((2008) . Such extremely low values of R e 
are visibly different from all other high-z data. They have 
median radii 2.5 to 3 times lower than t he global WINGS 
mass-r adius relation. It has to be noted that lvan Dokkum et alj 
(2008) extreme cases lay in a section of the plot where we are 
loosing our model resolution 8 in age. On the other hand, for 

8 It is very difficult to properly assign an age to galaxies older than 9Gyrs, 
as their spectra are practically the same. 



(2009). 

The present study demonstrates that comparing the sizes 
of passive galaxy samples at high-z with local samples, can 
mimic a fictitious evolution of radius with redshift, if the ef- 
fects of the stellar age selection are not properly taken into 
account. In constrast, our results point towards an overall con- 
sistency between the sizes of high-z quiescent and low-z old 
massive galaxies, with at most a very mild evolution in size, 
as far as cluster galaxies are concerned. Considering single 
points in Fig|4] instead of the global tendency, the maximum 
amount of evolution in size from the present study is a fac- 
tor of 1.5, much lower than the claimed factor 3. In a picture 
where systematic er rors affect high-z measures of sizes (see, 
Mancini et al. 2009) such a factor would be probably easily 
accomodated, but this will have to be further investigated. 
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FIG. 4. — Median R e of WINGS galaxies (filled dots) that stopped forming stars (i.e. with luminosity-weighted ages older by) at least > 1.5 Gyr before the 
redshift plotted. Three stellar mass intervals are considered. High-z literature data are plotted as color symbols (see legend). WINGS error bars represent the 
errors on the median size, while the dotted error bars are the upper and the lower quartiles of the corresponding distribution. Due to the scarcity of high-z data, 
the high-z datapoints are mean values with their corresponding RMS scatter (dotted error bars) and RMS of the mean (full error bars) When the mean is done on 
less than 3 galaxies, the symbol is dotted, to show its low statistical significance. 



such extreme cases systematic effe cts caused by large d i stance 
could be important; for example Ivan Dokkum et ail (120081) 
discuss some caveats on size estimates of their high-z sample 
which could give a factor o f 2 greater sizes, muc h more com- 
patible with FigH](see also, Mancini et al. 2009). However, it 
is true that, if more galaxies of such compact nature would be 
found in the future, they would be candidates of a "growing- 
radius" class of galaxies which would not be explained by an 
age selection effect. 

4.1. Frequency and number density 

We now turn to analyze the frequency and number density 
of WINGS SDGs. We have seen that SDGs represent a siz- 
able fraction (22%) of all cluster spectroscopically confirmed 
members more massive than 3xl0 10 M Q . This fraction does 
not vary using higher-mass cutoff limits, i.e. 5 or 8xlO ll) M . 

We determine the expected total number of SDGs in all 
WINGS clusters by multiplying the average SDG number per 
cluster among the 21 clusters considered in this study, cor- 
rected for spectroscopic completeness (~ lOSDGs/cluster), 
by the total number of clusters in the WINGS survey (78). 



We then calculate the whole comoving volume associated 
with the redshift range of WINGS clusters: 

Vwings = y (Rl -Ri) (1 - sinb) = 5.73xl0 7 Mpc 3 (5) 

where b = 20° is the limit in galactic latitude imposed by the 
survey to avoid the galactic disk regions, and Ri = 169.8 Mpc 
and R2 = 295.0 Mpc are the distances in our cosmology corre- 
sponding to the minimum (z = 0.04) and maximum (z — 0.07) 
redshifts of our clusters, respectively. 

Assuming no SDG is present outside of WINGS clusters 
in this volume, a very hard lower limit to the SDG number 
density in the local Universe is then N = 1.31xlO _5 Mpc~ 3 
for M, > 3xl0 10 M Q , and N = 046xlO- 5 Mpc~ 3 for M* > 
8xlO lo M (see TableHJ. 

Considering only the volume effectively probed by the 
WINGS clusters (a total area of about 25 sq.deg. and an av- 
erage redshift range z ± 0.007 around each cluster redshift), 
the SDG number density in clusters turns out to be very high, 
JV = 2.9xlO- 2 Mpc~ 3 . 

Ideally, we would like to compare the SDG number den- 
sity we derive with the number density of high-z SDGs, to 
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investigate what fraction of the distant SDGs can have super- 
dense local descendants, consistent with having maintained its 
size and mass unaltered since z > 1. Unfortunately, the SDG 
number density at high-z is not available. The information 
that several authors provide is the number density of high-z 
quiescent galaxies (ICimatti et al.f 2008: iBezanson et ai]l2009l: 
IWuvts et al]|2009l) . but, as it can be also seen in FigQ] a large 
fraction of high-z quiescent galaxies are not superdense. 

The number density of high-z quiescent galaxies can be 
compared with WINGS estimates for galaxies that according 
to our luminosity-weighted ages should be quiescent (= with 
stellar ages older than 1.5 Gyr) at each redshift (TableHJ. 

Interestingly, nearly 20% of all high-z quiescent galaxies 
are found in low-z WINGS clusters, of which about one third 
is superdense. Several uncertainties and systematics affect 
this fraction: for example, just by taking all galaxies older 
that 9 Gyr instead of 10 Gyr, the number density of quies- 
cent galaxies in clusters increases by 50%. Moreover, we 
stress that the number of both SDGs and quiescent galax- 
ies we find in WINGS are not a complete census in clusters 
at z = 0.04 — 0.07, for several reasons: a) WINGS typically 
probe out to about half of the cluster virial radii; b) WINGS 
is an X-ray flux-limited sample, with different flux limits in 
the Northern and Southern hemisphere, and it is not complete 
down to a fixed X-ray luminosity; c) WINGS does not include 
clusters with L x < 0.2xl0 44 erg/s. 

This is an indication that the remaining 80% of quiescent 
galaxies at high-z and, presumably, some fraction of the SDG 
population, must be found today in the "field", where the 
"field" includes the outer regions of WINGS clusters, as well 
as a large number of low-mass clusters and groups. Hence, we 
speculate that groups and clusters may host a large fraction or 
even the totality of high-z quiescent descendants. 

5. WHY LOW-Z SDGS GALAXIES WERE NOT FOUND BEFORE 

In Fig|5]we plot the mass-radius relation for M* > 1O 1O M 
WINGS cluster m embers (small black dots). The SDSS 
IShen e t al. (2003) relation, commonly used by high-z stud- 
ies as local reference, is superimposed as a red full line for 
early-type galaxies (Sersic n > 2.5), and as a blue dashed 
line for late-type (n < 2.5), with their ler limits. Large dots 
are WINGS median values with upper- and lower-quartiles 
(completeness corrected), blue open dots for late-type galax- 
ies (morphologies later than SO), and red filled dots for early- 
type (elliptical and SO) galaxies. Green squares are the Bright- 
est Cluster Galaxies (BCGs): together with all galaxies with 
M* > 4x10' 'Mq these appear to significantly deviate from 
the general trend of all galaxies, showing a marked steepen- 
ing of the mass-radius relation at high masses. Indeed they 
seem to belong to a separate class of objects which likely 
underwent a significant accreti on of cold gas during the for- 
mation of the cluster (see, e.g.. lBernardill2009l) : this was the 
main reason why we considered inappropriate to include them 
in the present study. We want to stress here that this is just 
an empirical upper mass limit that arises from a visual in- 
spection of Fig|5] whose physical explanation is simply ten- 
tative and qualitative. With this is mind, we also note that 
high-z data (symbols and color code are the same of Figj4j 
of M* > 4x10' 'M are smaller by a factor of > 3 than the 
WINGS cluster BCGs. A strong evolution in radius is thus 
required for BCGs, at odds with the rest of the galaxies. 

Importantly, we find systematically lower radii (~ 0.1 dex) 
in our cluster early- and late-type galaxies, when compared 
with SDSS. In Tabf2]we report, as reference, our median R e 



values for different bins in mass, both for early- and late-type 
galaxies, t ogether with the co rresponding SDSS value calcu- 
lated from IShen et ail (120031) . If low-z cluster galaxies are 
the proper descendants of the high-z ones, ~ 20 % of the 
claimed evolution in radius needed to match the local mass- 
radius relation could be due to the uncorrect choice of the 
local relation. It is worth noting that the difference in the lo- 
cal mass-radius relation could be almost completely explained 
with the systematic offset in mass we are finding w ith SDSS- 
DR7 mass es, at variance with [ Gallazzi et al. (2005), Hyperz- 
mass and lBell & de Jongl(l2001l) (see section 13. 11 1. 

In the following we briefly discuss the possible causes of 
systematic errors which one might take into consideration 
when assessing the problem of how and how much high-z 
galaxies have undergone structural evolution. 

Galaxy stellar masses. As discussed in §2, when compar- 
ing different datasets it is of paramount importance to ensure 
consistency on the IMF assumed and, at high-z, on the model 
prescriptions. IMF slopes and limits have to be carefully 
matched. At high-z, it is also important to homogenize the 
treatment of the stellar TP-AGB phase, since the masses of 
stellar populations with ages of approximately 2Gyrs, can be 
over- estimated by ~ 0.15 dex (Maraston 2005; Cima tti et alj 
2008). Furthermore, the type of stellar mass considered is 
crucial (see §2). At high-z, the mass locked into stellar rem- 
nants is negligible, while it becomes more and more important 
at lower-z (up to ~ 0.15 dex). Hence, the comparison with 
low-z masses, to be meaningful, should be carried out consid- 
ering the mass locked in remnants at low-z (masses n.2 in §2). 

Effective-radius. We note that iBlanton etaTJ (120051) dis- 
cuss a bug in the 2003 tool which measured the sizes for 
the NYU-VAGC catalog, that caus ed small radii to b e over- 
estimated; these sizes were used bv lShen et al.l d2003l) for the 
SDSS mass-radius relation. While we are not able to estimate 
the importance of this effect, the corrected radii should clearly 
be used to reassess the SDSS mass-rad ius relation. Further 
more, a recent paper bv lGuo et al.1 (|2009), studies the possible 
biases induced by a no isy background subtra ction. 

In a recent paper, iMancini et al.1 (120091) claim that R e 
measurements of low S/N high-z compact galaxies, may give 
systematic lower sizes up to a factor of ~2. 

Extreme selection criteria. The definition of a "su- 
perdense galaxy" i s necessarily arbitrary at some level. 
iTrujillo et al.1 J2009h search for SDGs with M* > 8xl0 10 M Q 
and R e < 1.5kpc, and found no candidates with old 
luminosity-weighted ages, and very few of all ages. Only 
9% of the high-z galaxies (10 out of 108) considered in our 
analysis fulfill these extreme mass and radius selection crite- 
ria. Interestingly, we found approximately the same fraction 
(~ 5%), applying this definition in our sample: we have only 
16 "extreme" SDGs, with a median luminosity- weighted age 
of 10. 1 Gyr. Hence, it is plausible that only a small fraction of 
galaxies satisfies such extreme criteria. 

Completeness. In a recent paper, iTaylor et ail ([2009) 
perform a thoroughly search of SDGs in the z ~ 0.1 SDSS 
galaxies. In particular the authors discuss the different 
varieties of incompletenesses involved when using SDSS 
data to assess a proper local mass-radius relation. Fiber 
collision and/or limit in surface brightness when selecting 
galaxies for spectroscopic follow-up, could cause SDSS to 
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TABLE 1 

Number densities o f SDG s and q uiescent galaxies. L iterature data: Bez= IBezanson et al.I j2009h . Cimatti= ICimatti et al.I 
j2008l) .WuYTs MWuYTS et alU 2009). Errors are derived prom Poissonian statistics. 



Criteria 



WINGS Literature 

10~ 5 Mpc~ 3 10~ 5 Mpc~ 3 



SDGs 

SDGs M* > 8xl0 10 M Q 
SDGs Lw - age > lOGyr (z = 1 .5) 
SDGs R c < 1.5kpc 
Quiescent, 3xl0 10 M Q < M* < 4xlO u M 

Quiescent z ~ 2.5, M* > 10" M Q 
Quiescent z ~ 1.5, 10 10 < M* < 10"Mq 
Quiescent z > 1.5, M* > 4xlO'°M 
Quiescent z > 1.5, M» > 10" M 



1.31 ±0.09 
0.46 ±0.05 
0.57 ±0.06 
0.68 ± 0.07 
1.55 ±0.06 - 
0.50 ±0.06 
1.66 ±0.10 
1.80±0.11 
1.09 ±0.08 



Bez=3.6 
Cimatti=10 
Wuyts=ll 
Wuyts=4.5 



TABLE 2 

Local WINGS mass-radius relation. The values are 

LOGARITHM OF THE MEDIAN ESTIMATES, ERRORS ARE THE LOWER AND 
UPPER QUARTILES OF THE DISTRIBU TIONS. THE SPSS Re IS THE 
EXPECTED VALUE CALCULATED FROM lSHEN ET AL.I j2003) AT THE 
SAME MASS REPORTED IN THE FIRST COLUMN. 

log 10 (M./M B ) log 10 (R e /kpc) log 10 (R e SDSS /kpc) 
Late Type galaxies 
9^F^ 0lF^ 041 
10.27+°^ 0.40t°-;j 0.48 
10.571°^ 0.47+°^ 0.55 

10.86too7 - 48 to'i3 063 

11.06+.°,;!® 0.67 '+!■'" 0.70 



Early Type galaxies 
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miss local clusters galaxies in a systematic way. However, 
the authors show that completeness issues are not sufficient 
to explain the lack of ga laxies as compact as those in the 
Ivan Dokkum et al] ([2008!) sample, which do not have a local 
counterpart in WINGS too (see section HJ. We refer to that 
paper for the details on their analysis of massive compact 
galaxies in the SDSS. 

6. CLUSTER SDGS PROPERTIES: CLUES TO THEIR ORIGIN 

The WINGS SDGs sample consists of 31 ellipticals, 78 
SOs and 13 late-type galaxies. When completeness corrected, 
these numbers become 46.5 (22.8 ±4%), 136.9(67.3 ± 7%) 
and 20. 1 (9.9 ± 2%), respectively. If compared with the overall 
morphological fractio ns in a magnitude limited sample (see 
Poggi anti et 10120091 Fig.l and Tab.l), there is an excess of 
SOs at the cost of ellipticals and later types (see Tab[3]). This 
might indicate either that the SO morphology is preferred by 
SDGs, or that some of these SOs have uncorrect R e because 
their light profile is not well suited for a single Sersic law fit. 

In Fig|6] we present the distributions of several relevant 



TABLE 3 

Characteristic numbers of WINGS SDGs. Errors are derived 

FROM POISSONIA STATISTICS FOR COUNTS, AND ARE RMS FOR OTHER 
QUANTITIES. C.C. COMPLETENESS CORRECTED. 



Quantity 


Value 


RMS eiror 


SDGs 


122 


11 


SDGs C.C. 


203.5 


14.3 


<Re> 


1.61 


0.29 


{„) 


3.0 


0.6 


(b/a) 


0.54 


0.18 


<M*> 


8.7xlO lo M 


2.5xlO'°M 


(Vabs) 


-20.68 


0.37 


(Lw - age) 


9.62 


1.94 


(Mw — age) 


12.02 


1.28 


Ellipticals frac. C.C. 


22.8% 




SOs frac. C.C. 


67.3% 




Late-type frac. C.C. 


9.9% 





quantities describing our SDGs sample. First of all, the axial 
ratio distribution ({b/a) = 0.54 ±0.18) shows that WINGS 
SDGs have a tendency to be flattened, mostly due to the high 
fraction of SOs. As expected, the population of elliptical 
galaxies is remarkably rounded. The late-type galaxies are 
extremely flat and could, in principle, introduce systematics 
in our analysis. We decided to keep them because most of the 
high-z samples are not selected on the basis of their morphol- 
ogy and therefore may include late-type galaxies. Anyway, 
all the conclusions of our present study are even reinforced if 
only early-type galaxies are considered. 

The Sersic indexes of our SDGs are characteristic of disky- 
like rather than of early-type galaxies ((«) = 3 .0 ± 0.6): while 
late-type SDGs present an expected n ~ 1 value, the majority 
of elliptical SDGs have remarkably low values too. 

The WINGS median R e is similar to that of high-z SDGs 
((R e ) = 1.61 ±0.29), with a few of them being larger than 
2.5 kpc. These objects are the most massive ones, and they 
probably are transition objects from the compact phase to a 
more complex radius inflation phase, where most probably 
galaxies acquire gas and/or stars in the external regions, in- 
creasing R e (see Fig|5]l. 

Our SDGs have high intrinsic luminosity (My) = —20.68 ± 
0.37, and stellar masses (M*) = (8.67 ± 2.55)xl0 10 M Q . 

WINGS SDGs may show a slight tendency to prefer the 
central regions of clusters (CCD) and intermediate/high den- 
sity regions (LD), but overall their clustercentric and local 
density distribution are not too dissimilar from those of galax- 
ies of similar mass. Our images are uniformly sampling the 
cluster populations as far as R^q (~ 0.6,/? 200) f° r a U clusters, 
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FIG. 5. — Mass-radius relation for WINGS cluster members (small black dots). Open blue circles are the median values with upper- and lower-quartiles for 
late-type (later than SO) galaxies, and red filled circles for e arly-type (elliptical and SO) galaxies. Green squares are WINGS Brightest Cluster Galaxies (BCGs). 
The SDSS median and lcr relations from Shen et al. (2003) for early- (red full line, n > 2.5) and late-type (blue dashed line, n < 2.5) galaxies are also drawn. 



so the sharp decline at larger radii may be just a result of the 
area coverage. A future ancillary project with the forthcom- 
ing OMEGACAM at the VST telescope will survey a consid- 
erable fraction of WINGS clusters at much larger radii, and 
will uncover possible compact candidates at larger distances. 

We have already discussed the ages of our SDGs in §3, 
showing the difference between Lw-age and Mw-age, and 
the care that needs to be taken when using these quantities 
to select samples. WINGS SDGs have high median Lw- 
((Lw — age) = 9.6 ± 1.9Gyr) and Mw-ages ((Mw-age) 
12.0 ■ 1.3Gyr), showing that the bulk of the mass has an age 
typically 2 Gyr older than the luminosity-weighted age (red 
dashed histogram in Fig|6]l. 

To quantify to what extent, on average, WINGS SDGs are 
older than non-compact galaxies of similar masses, we used 



the Monte Carlo technique to extract 1000 random samples 
of "normal" galaxies with the same mass distribution of the 
SDGs. This is done to disentangle the dependence of age from 
mass. We used the early-type galaxies only, to be more con- 
servative, as we know that there are more late types among 
"normal" galaxies than in the SDGs sample. In Fig|7]we plot 
the distributions of the median Lw- (blue lines) and Mw-ages 
(red lines) of these 1000 samples. Choosing "normal galax- 
ies" according to the SDG mass distribution gives on aver- 
age younger ages than those of all (i.e. not mass matched) 
"normal" galaxies (vertical dashed lines). Importantly, the 
Monte Carlo simulation shows that compact galaxies tend to 
be ~1.5Gyr older (both in Lw- and Mw-) than normal galax- 
ies of the same mass, again suggesting that in some way age is 
related to compactness, in addition than mass, as we discussed 
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FIG. 6. — Distributions of different quantities of interest for WI NGS SDGs clus ter members. The color shades correspond to morpholo gies, red=ellipticals , 
green=SOs and blue=late-type galaxi es. M t = total stellar mass I Fritz et al. 2007), CCD is the cluster centric distance in units of R200 (Carlberg et al J l 19971) . 
while LD is the local density I Dressier 1980). The magenta dashed histograms are the distribution of non-compact galaxies in the same mass range of SDGs. 



regarding Fig. 2. 

We have seen that the WINGS SDGs sample consists of 
galaxies similar in all respects to the compact quiescent ones 
found in recent high-z studies. They are very old, massive, 
and compact. Their presence in the local universe opens 
new perspectives on their formation and evolution, which may 
change the present understanding of their nature. 

7. SUMMARY AND CONCLUSIONS 

We find 122 SDGs in the WINGS survey of nearby galaxy 
clusters (z ~ 0.05), with stellar mass 3xl0 10 < M*/M© < 
4xlO n and surface mass density £50 > 3xl0 9 M©kpc -2 . They 
represent nearly 22% of all cluster members in the same mass 
range. They have masses and sizes similar to their high-z 
counterparts. 

We find that both mass and radius determine the age of 
massive low-z cluster galaxies: the larger the mass, and the 
smaller the radius, the older the stellar population. Select- 



ing quiescent galaxies at any redshift results in selecting the 
smallest galaxies; the further back in time we search for qui- 
escent galaxies, the smaller the sizes we measure as a conse- 
quence of this effect. 

We compare our data with spectroscopic high-z studies, 
whose mass and size estimates are more reliable than pho- 
tometric ones. If cluster galaxy sizes and masses today are 
representative of those of high-z galaxies, our findings show 
that there is no need for an evolution in size (at least as far 
as z ~ 2), when this age effect is properly taken into account. 
The largest possible evolution in size at z < 2, based on dy- 
namical and therefore possible dishomogeneous mass mea- 
surements, is very mild, a factor of 1.5 at most, much smaller 
of the factor 3 claimed in the literature; anyway, it is difficult 
to directly interpret this discrepancy due to the different way 
masses are estimated. In contrast, there is strong evidence for 
a large evolution in radius for the the most massive galaxies 
with M* > 4x10' M© compared to similarly massive galax- 



Superdense galaxies in WINGS clusters 



13 



■ Lw-Ages 




6 8 10 12 



Age (Gyr) 

FIG. 7. — Distributions of median luminosity-weighted ages (left blue his- 
togram) and mass- weighted ages (right red histogram) of 1000 random ex- 
tractions of non-superdense early-type galaxies with the same mass distribu- 
tion of SDGs. The vertical lines are the median ages of SDGs (solid lines) 
and of non-superdense galaxies (dashed lines) with their intrinsic mass dis- 
tribution. 

ies, in WINGS, i.e. the BCGs. 

On the other hand, the sizes of galaxies in the sample of 
Ivan Dokkum et al] (|2008) at z = 2.4 are smaller by a factor of 
~ 3 even with respect to the WINGS datapoints in FigJU when 
the age-selection effect is taken into account. These objects 
represent a population of galaxies visibly different from other 
high-z data (see FigQ}. 

For masses M„ < 3xl0 I0 M Q there is a significant decrease 
in the frequency of SDGs and speculate that this rapid de- 
crease in number could be an indication that a minimum 
threshold in mass is required to form compact galaxies . 

The local mass-radius relation by IShen et al.l d2003i) . used 
by high-z studies as reference, turns out to be shifted toward 



higher radii at fixed mass when compared to the WINGS rela- 
tion. This is probably due to the systematic offset between our 
total masses with respect to SDSS-DR7 masses, discussed in 
section [3~T1 anyway, o ur masses turn out to be in good agree- 
ment with SDSS-DR4lGallazzi et all (120051) . Hyperzmass and 
iBell&de Jong (120011) . 

Assuming that SDGs reside only in clusters, we calculate a 
hard lower limit of their number density in the nearby universe 
of 1 .3xl0~ 5 Mpc~ 3 , which becomes 0.57xl0~ 5 Mpc~ 3 if only 
SDGs that were quiescent atz ~ 1.5 (i.e., luminosity- weighted 
age > lOGyr) are considered. While no published data on 
high-z SDGs number density is available, there are estimates 
for the density of quiescent galaxies at z ~ 1.5. We find a 
lower limit of 0.18xl0~ 4 Mpc~ 3 of such quiescent galaxies 
with M* > 4xlO 1() M0 in clusters, to be compared with a high- 
z value of 10 _4 Mpc -3 (IWuyts et alJl2009h . Around 20% of 
all high-z quiescent galaxies are therefore found in the inner 
regions of WINGS clusters. 

Our findings challenge the simple picture of a widespread 
evolution of the radius of compact high-z galaxies with red- 
shift. The presence of compact galaxies in local clusters sug- 
gests that the formation and evolution of such systems may 
not be simply explained with a "growing radius" mechanism, 
as presently thought. In particular our research can be used to 
further constrain the current picture of galaxy mass assembly 
(hyerarchical merging, down-sizing, etc.) from the first few 
Gyrs after the Big Bang, to the present galaxy clusters and 
high density regions probed by the WINGS survey. 
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